Time dynamic of Fourier modes in turbulence: 
Sweeping effect, long-time correlations and time intermittency 
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We present an experimental study of the typical time scales of spatial Fourier modes of the 
vorticity in turbulence. The measurements of time auto-correlation functions, for well defined spatial 
wavenumbers k, reveal that vorticity exhibits two characteristic times: a short one depending on 
the inverse of the wavenumber, which can be interpreted as a manifestation of the sweeping effect 
(random advection), and a larger one, independent of k, of the order of the integral time scale of the 
flow. In a second step, we show that two different vorticity Fourier modes with wavenumbers k and k' 
are long-time cross-correlated, whatever the spectral gap k — k'. The existence of a significant cross- 
correlation between scales, indicates a strong statistical dependance between scales, at variance with 
the original Kolmogorov theory (1941). This long-time correlation, closely related to the integral 
forcing time, could be ascribed to some evidence of time intermittency, associated with fluctuations 
in time of the energy injected at large scales and resulting in a lack of global scale invariance of the 
dynamical properties of the turbulent flow. 

PACS numbers: 47.27.Gs, 47.32. Cc, 43.28.-g 



Statistical intermittency, responsible for a lack of sta- 
tistical global scale invariance, remains one of the main 
puzzling features of turbulence, still unresolved Q. At 
the experimental level, following Kolmogorov initial pre- 
diction 0, the evolution across scales of the turbulence 
statistical properties can be traced by computing the spa- 
tial velocity increments, 5u r (x,t) — u(x + r,t) — u(x,t). 
Usually, one resorts to hot-wire anemometry providing 
Eulcrian measurements of the longitudinal velocity com- 
ponent at one point x along time t. Scale dependence 
is then recovered by mapping time increments St of the 
Eulerian signal onto spatial increments r = —U avg 8t ac- 
cording to the Taylor hypothesis of frozen turbulence : 
5u r (t) = u{x 0l i) — u{x ,t — r/U a vg)- In this context (spa- 
tial intermittency), and according to Kolmogorov refined 
similarity hypothesis statistical intermittency is the 

consequence of the strong spatial heterogeneity (possibly 
multifractality) of the local energy transfer rate, related 
to some kind of multiplicative cascade process of the tur- 
bulent energy across scales. In the spirit of Kolmogorov 
approach, the energy transfer rate is also expected to dis- 
play strong time fluctuations, reflecting unavoidable fluc- 
tuations of the energy injected at large scales. Recently, 
it has been suggested that the study of dynamical 
multiscaling (time dependence of the velocity structure 
functions) could be of valuable help for a better compre- 
hension of turbulence intermittency. Unfortunately, as 
implicitely stated by the Taylor hypothesis, the statistics 
of Eulcrian velocity increments are strongly dominated 
by the spatial features of the flow and thus, weakly sen- 
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sitive to its time fluctuations. Actually, it is generally 
accepted that the latter fact is a direct consequence of 
the random character of the advection by the large scale 
flow, past any Eulerian probe, of the whole velocity field 
(the so-called sweeping effect [H|(J). 

An alternate and efficient way to study statistical scale 
dependence is to rely on a spatial Fourier analysis of the 
field of interest UQ- Spectral measurements, based on 
wave scattering experiments are commonly used in var- 
ious domains of research in physics, starting with phase 
transitions and critical phenomena in condensed matter 
physics where light and neutron scattering techniques 
are largely widespread. It is now well established on 
both theoretical j^, an d experimental basis 1 that 
acoustic waves propagating in a turbulent medium can 
be scattered by vorticity fluctuations. As in any scatter- 
ing experiments, it may be shown 0] that, due to the 
coherent average of the waves scattered by the vorticity 
distribution, the overall scattered amplitude p sca tt(t) is 
linearly related to the spatial Fourier transform of the 
vorticity field according to : 



Pscattit) CC LO±_ (k,tj .po(t) 
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the scattering wavevector k being a function of the 
incoming sound frequency vq and the scattering angle 9: 



k = Aiyisq sin (9/2) jcx 



(3) 
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FIG. 1: Sketch of the experimental setup: turbulence is pro- 
duced by a turbulent round jet at R\ ~ 600. The acoustic 
measurement volume V SC att is located downstream the noz- 
zle. The wavevector k is parallel to the mean flow axis (x) 
while the probed component of the vorticity vector field is 
perpendicular to the scattering plane. 
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FIG. 2: Time auto-correlation function Cfc(r) at wavenumber 
k — 1200m -1 . The global shape is similar for all analysed 
wavenumbers: a short-time decorrelation with characteristic 
time ts (see insert for detail), followed by a much slower de- 
crease, with (see also Fig. 0] (a) for a log-lin plot) and char- 
acteristic time tl- Here, ts = 0.88ms and tl = 45ms. 



with c the sound speed. Note that only the component of 
the vorticity vector field, normal to the scattering plane, 
is involved in the scattered amplitude (hence the index 
The scattering setup consists in a bistatic configura- 
tion (Fig. ^): an ultrasound monochromatic plane wave 
Po{t), with frequency vq, continuously insonifies the tur- 
bulent media and the acoustic amplitude p sca tt(t), scat- 
tered in the direction 9, is recorded along time by a re- 
ceiver. Both acoustic transmitter and receiver work in 
a linear regime (they are phase sensitive). According to 
equation l[TJl. a direct image of the spatial Fourier mode 
of the vorticity at wavevector k, is obtained by a sim- 
ple heterodyne demodulation providing a complex signal 
(phase and amplitude). 

The main interest of the acoustic scattering experiment 
lies in the fact that, contrary to classical Eulerian mea- 
surements, the scale dependence (through the selection 
of a unique wavevector k) is performed independently of 
the time evolution of the collected signal, thanks to the 
direct spatial Fourier transform involved in the scatter- 
ing process. We will now focus on the time behavior of 

the Fourier modes of vorticity w {k, tj , as a function of 
the length scale parameter k. Although a Fourier mode 
is a complex quantity (as is the demodulated scattered 
signal), we will restrict ourselves to the amplitude of the 

signal: (k, tj | hereafter noted u>(k, t) for sake of sim- 
plicity. We have investigated the statistical properties of 
a turbulent round air jet at room temperature. The di- 
rection of the probed vorticity field is radial, while the 
direction of the scattering wavevector k is aligned with 
the mean flow velocity. The flow emerges in the x direc- 
tion from a nozzle of diameter D = 0.12m (see Fig. QJ. 
Throughout the experiment, the scattering angle is kept 
at a constant value and different wavevectors are anal- 



ysed by varying the incoming sound frequency v . As 
the scattering angle 9 is constant, one can show that, in 
our bistatic configuration, the spectral resolution is given 
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by 5k ~ V S catt > independent of the analysed wavenum- 
ber k fl2i]. The measurement volume V sca tt is defined 
by the intersection of the incident and detected acous- 
tic beams and mainly depends on 9 and on the size of 
both acoustic transducers. With 9 = 40° and a diam- 
eter of the circular transducers of 14cm, the linear ex- 
tension of V scat t is of the order of the integral length 
scale of the jet flow estimated to L = 0.36m. Additional 
flow parameters have also been estimated, using conven- 
tional hot-wire anemometry: the Taylor micro-scale is 
At = 7.6mm and the associated Taylor-based Reynolds 
number R\ = u'\t/v worths R\ ~ 600 where v! is the 
longitudinal root-mean-square of the longitudinal veloc- 
ity fluctuations. 

From the time signals oj(k,t), collected at a fixed 
wavevector k, we compute the time auto-corrclation func- 
tion Ck(r) =< uj(k, t) u(k, t — r) > t where < . > t stands 
for the time average. A typical evolution of the normal- 
ized Cfc(r) (Cfc(0) = 1)), with respect to the time lag r, is 
sketched on figureEl Whatever the turbulent scale k, the 
same global shape is found, exhibiting two different and 
well separated characteristic times. For time lags close 
to zero, one observes a rapid decrease, with a more or 
less Gaussian shape. At larger time lags, a much slower 
decrease is visible, with a nearly exponential behavior 
as evidenced on the semi-log representation in Figure 01 
(a). Let us estimate the short time ts, by measuring the 
half amplitude width (Cfc(rs) = 1/2 in the small lags re- 
gion, and the long time t l , with an exponential fit of the 
tails of Cfc(r) at large times. Although such a Gaussian 
shape at short times has been predicted in some theo- 
retical models 0, ^| , the long time behaviour does not 
seem to have ever been reported. 
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FIG. 3: Loglog plot of the evolution of the two characteristic 
times of Fourier modes in turbulence. 



We now turn to the scaling of those two characteristic 
times with the wavenumber k. By tuning the incoming 
sound frequency (10kHz <v a < 170kHz), we have suc- 
cessively probed a decade of wavenumbers k, spanning 
the inertial range of the turbulent flow down to the dis- 
sipative range (0.30 < /cAt/27t < 4.7). Figure |3] clearly 
reveals two different behaviours: t$ is scale dependent 
while tl is not. Systematic studies, performed in var- 
ious flow configurations (velocity v! and integral scale 
V) gave us reliable results for both times. As for the 
short time ts, it decreases with k following the power 
law Tg oc (fcit')" 1 . In our experiments, the proportional- 
ity constant is about 1, and slightly decreases with the 
Reynolds number. Such a scaling law, involving the root 
mean square velocity (a large scale quantity), is usually 
ascribed to a sweeping effect associated with the random 
advection of the vorticity field by the large scale velocity 
fluctuations, according to Navier-Stokes equation. The 
consequence of the sweeping, with respect to Eulerian 
turbulence statistics, is a controversial question as re- 
gards the —5/3 exponent of thepower-law scaling of the 
turbulent energy-spectrum El El • However, 
it is generally well accepted that the sweeping effect is 
responsible for the unexpected (according to the global 
scale invariance hypothesis) f~ 5 ^ 3 power law behaviors 
of the spectra of powers (p > 2) of the Eulerian tur- 
bulent velocity fluctuations \u'\ p (where / is the tempo- 
ral frequency) [f| . Accordingly, our experimental results 
confirm that this latter effect can be more clearly evi- 
denced in the Fourier domain as suggested in reference 
[l8j |. Let us now turn to the large time tl- as we have 
checked extensively (by varying acoustic as well as hy- 
drodynamic conditions separately) tl is not scale depen- 
dent but rather depends on the forcing scales. Actually, 
we found that: tj, oc (/clw') _1 where k^ = 2ir/L is the 
large scale wavenumber. The fact that this long time 
correlation is not scale dependent and, besides, remains 
almost identical whatever the probed scale, with a sig- 



nificant correlation level (about 20% Figs. 0] (a)-(b)) is 
noteworthy. First, the time tl is exactly proportional to 
the integral Eulerian time (extracted from the velocity 
auto-correlation fonction see Fig. 01(c)), with a propor- 
tionality constant of one: tl w (fcz,u') . Second, the 
long-time invariant behaviour through the scales, con- 
trary to the scale dependence of the short time, indicates 
that it is not possible to find a scaling transformation 
which leads to a collapse of all the correlation functions 
onto a single one. We are lead to the conclusion that the 
coexistence of these two characteristic times indicates a 
lack of global scale invariance in the cascade process. In 
the spirit of spatial intermittency, related to the absence 
of an universal shape of the probability density functions 
of the spatial velocity increments, it is thus tempting to 
term our observations as a manifestation of time inter- 
mittency. It is worth mentionning, at this point, that 
we have also evidenced some influence of the spatial in- 
termittency of turbulence, manifested by a significant de- 
pendence of the statistics of the spatial Fourier modes on 
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the spectral resolution 8k or the spatial extension V sca j. t 
of the experimental setup with respect to the integral 
length scale of the turbulent flow jlflj . 

The close resemblance in the long time evolution of any 
pair of spatial Fourier modes with different wavenumbers 
k and k', suggests that they could be dynamically driven 
by the very same large scale process. To get further in- 
sight about this observation, we have performed simulta- 
neous acquisitions of two distinct spatial Fourier modes k 
and k'. Provided the spectral gap k — k' is large enough, 
such measurements can be performed by driving a single 
transmitter with the sum of two sine waves with the ap- 
propriate frequencies v a and v' Q (Eqn. |2J. The scattered 
pressure signals, around each incoming frequency, can be 
easily separated by means of two simple band-pass fil- 
tering operations. To avoid spurious interference effects, 
the same investigation could also be performed by using 
a second pair of transducers defining a second indepen- 
dent scattering channel as in |2jj. However, the single 
pair configuration presents the advantages of a better 
wavevector alignment as well as the best possible mea- 
surement volume matching. Actually, we have carefully 
checked that both setups give the same results. From two 
synchronous time series (at k and k ), with a well con- 
troled spectral gap k — k' ', we have computed the cross- 
correlation function C^wir) =< u>(k, t) u>(k', t — r) > t , 
of two spatial Fourier modes. A typical example of this 
cross-correlation is represented on Fig. 0] (c) (the same 
behaviour is observed whatever the scale separation) . We 
have also represented on Fig. Q](a) (resp (b)), the auto- 
correlation function of the two spatial Fourier modes with 
wavenumbers k (resp. k'). A significant level of cross- 
correlation is observed for time lags up to the integral 
time scale, with a shape and an amplitude close to that 
of both auto-correlation functions. The main difference 
between auto-correlation and cross-correlation functions 
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FIG. 4: Log-lin plots of time correlation functions: (a) auto- 
correlation Cfe(r) for k = 1200m -1 (as in Fig. |2J, (b) auto- 
correlation C£.(t) for k — 820m -1 , (c) solid line: cross- 
correlation Cfefc'(r) between the two foregoing Fourier modes, 
dotted line: auto-correlation C uu (j) of the Eulerian longitu- 
dinal velocity u at the center of V 3C att- 



at smaller Reynolds numbers (in a grid turbulence at 
Re\ ~ 100), as well as at much higher Reynolds numbers 
(up to Re\ ~ 6000, in a cryogenic Helium jet facility at 
CERN [221]), this could well be a generic feature of flow 
turbulence. 
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lies in the absence of the rapid decay at small time lags for 
the cross-correlation. We are lead to the conclusion that 
the two observed behaviours (short and long time corre- 
lation) should be associated to two different dynamical 
process driving the vorticity field. We also point out the 
fact that the existence of a significant cross-correlation 
between the two spatial Fourier modes k and k' implies 
a strong statistical dependance between scales, irrespec- 
tive of the spectral gap k — k'. This new result, is clearly 
at variance with the initial Kolmogorov 1941 theory [2J. 
As far as the long time driving process is concerned, the 
significant cross-correlation level suggests that all scales 
are instantaneously driven by the same process (possibly 
multiplicative) . Our belief is that the observed long time 
statistical dependance could be related to some memory 
effect of the time fluctuations of the energy injection at 
large scales, in the spirit of Landau's objection formu- 
lated against the K41 model in 1944 [lll2l|. 

To summarize, acoustic scattering allows the direct 
spectral probing, continuously in time, of spatial Fourier 
modes of one component of the vorticity vector field. 
Thanks to a good spectral resolution, the proper selection 
of a well defined wavevector k of the turbulent flow results 
in an unambiguous separation of the spatial and time 
features of the turbulent dynamic. As a first result, we 
clearly evidence sweeping effects, related to the random- 
ness of the large scale velocity field advecting the small 
scales and manifesting as a rapid decorrelation with a 
nearly Gaussian shape (reflecting the Gaussian statistics 
of the spatio-temporal advecting velocity field). A much 
slower modal decorrelation, over time lags up to the inte- 
gral time scale, is also observed and reveals a significant 
statistical dependency between scales. We propose to as- 
cribe our experimental observations to some large scale 
time intermittency (in contrast with the usual small scale 
spatial intermittency affecting the statistics of Eulerian 
spatial velocity increments). A plausible explanation for 
the source of such time intermittency, could be the ran- 
dom nature of the time fluctuations of the rate of energy 
injected at large scales. To conclude, we want to stress 
that, as we have observed similar statistical behaviours 
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